Quench dynamics of the topological quantum phase transition in the Wen-plaquette 
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We study the quench dynamics of the topological quantum phase transition in the two-dimensional 
transverse Wen-plaquette model, which has a phase transition from a Zi topologically ordered to a 
spin-polarized state. By mapping the Wen-plaquette model onto a one-dimensional quantum Ising 
model, we calculate the expectation value of the plaquette operator Fi during a slowly quenching 
process from a topologically ordered state. A logarithmic scaling law of quench dynamics near the 
quantum phase transition is found, which is analogous to the well-known static critical behavior of 
the specific heat in the one-dimensional quantum Ising model. 



I. INTRODUCTION 

Ultracold atoms provide an ideal platform for exper- 
imental studies of the time evolution of quantum sys- 
tems, and make it desirable for related theoretical explo- 
rations on dynamics of quantum phase transitions in var- 
ious models. These explorations mainly focus on nonequi- 
librium dynamics in quantum systems which undergo a 
quantum phase transition when a system parameter is 
varied (quantum quench) 1 . These experimental and the- 
oretical studies can potentially help to pave the way for 
future technologies and provide a deeper understanding 
of quantum many-body physics, particularly the univer- 
sal scaling behavior in the quench dynamics. 

Recently, a new type of phase transition, the so-called 
topological quantum phase transition (TQPT) has at- 
tracted considerable research attention 2-13 . TQPT is a 
kind of phase transition between two quantum states 
with the same symmetry. It is fundamentally different 
from the usual symmetry-breaking phase transition, and 
involves a new type of order — topological order, as in- 
troduced by Wen 14 . In such an ordered state, there is no 
local order parameter, and the state is robust against 
arbitrary local perturbations. On this basis, quantum 
systems with topological order have been proposed to 
build robust quantum memories 6 and topological quan- 
tum computer (TQC) 15 ' 16 . In Refs. 15,16 , it was shown 
that the pure state of topological order can be obtained 
via an adiabatical and continuous evolution process from 
a non-topologically ordered state and further be used as 
the initial state for TQC. Nevertheless, the detailed dy- 
namics of such a quench process, particularly the univer- 
sal scaling behavior near TQPT, has not been studied 
yet. 

In the last decade, several exactly solvable spin mod- 
els with topological order were found, such as the toric- 
code model 6 , the Wen-plaquette model 17 and the Kitaev 
model on a hexagonal lattice 18 . These spin models pro- 
vide a framework to study the TQPT and its quench 
dynamics. Thus, in recent years, some research groups 
have studied the quench dynamics of TQPT in the Kitaev 
model (e.g. Mondal et al. 19 ) or the toric-code model (e.g. 



Tsomokos et al. 20 ). Mondal et al. found a relationship 
between the quench rate and the defect density in the 
one-dimensional (ID) and two-dimensional (2D) Kitaev 
model in the limit of slow quench rate and generalized the 
result to the defect density of a d dimensional quantum 
model. Tsomokos et al. investigated how a topologically 
ordered ground state in the toric-code model reacts to 
rapid quenches. They tested several cases and showed 
which kind of quench can preserve or suppress the topo- 
logical order. 

In this work, we study the quench dynamics of TQPT 
from a topologically ordered state to a non-topological 
order in the Wen-plaquette model. To characterize the 
phase transition, we calculate the expectation value of 
plaquette operator F, during the quenching process, 
which is related to the number of quasiparticles in the 
topologically ordered state. Our results provide helpful 
information about the whole quenching process. 

The remaining of this paper is organized as follows. 
Section II describes an exact mapping from the 2D trans- 
verse Wen-plaquette model onto the ID Ising chain. In 
Sec. Ill, we study the TQPT of the 2D transverse Wen- 
plaquette model and give some results about its order 
parameters. In Sec. IV, we give the solutions to the dy- 
namics of TQPT in the transverse Wen-plaquette model. 
A brief discussion is given in Sec. V. 



II. MAPPING THE TRANSVERSE 
WEN-PLAQUETTE ONTO ISING MODEL 

We start with the Hamiltonian of the Wen-plaquette 
model on a square lattice with periodic boundary condi- 
tions in both directions: 
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(1) 



where rf and rf are Pauli operators on site i, x and y 
are the unit vectors in x-axis and y-axis, respectively (see 
Fig. 1). Because of the commutativity of H and Fi, the 
energy eigenstates can be labeled by the eigenstates of Fi. 
We can easily find Ff = 1, so the eigenvalues of Fi are 
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FIG. 1: (Color online) The lattice where Wen-plaquette model 
is located. A plaquette is defined by Fi = rf Wf +e .T? +&+ *Tf + $. 

i 71 ; = 1 and Fi — — 1, which gives the exact ground state 
energy. In the case of g > 0, the ground state is Fi = 
1 for every plaquette, and the elementary excitation is 

= -1 on one plaquette (denoted by i) with an energy 
gap Eg — Eq = 2g. On an even-by-even lattice, there are 
two types of plaquettes — the even plaquettes and the 
odd plaquettes respectively. As a result, one may define 
two kinds of bosonic quasiparticlcs: Z2 charge and Z2 
vortex (see detailed calculations in Ref. 17 or Ref. 2 ). A Z 2 
charge is defined by Fi = — 1 on an even sub-plaqucttc 
while a Z2 vortex defined by Fi = — 1 on an odd. Thus 
a fermion can be regarded as the bound state of a Z 2 
charge and a Z2 vortex. 

Now we consider the Wen-plaquette model in a trans- 
verse field, which is defined by: 

H' w = -gJ2Fi-jJ2rr- (2) 

i i 

This model on a square lattice can be mapped onto the 
ID quantum Ising model with the Hamiltonian 11 

jv 

Hj = -J2(9i<+«), (3) 

n=l 

where a„ and er* are Pauli operators. 

To derive the mapping, one can calculate the commu- 
tation relations (see detailed calculations in Appendix): 

[Fi , Tj ] — *2FiTj (Si^j — x ~1~ y)) 

[Fi,Fi\=0, 

[Tf,T/]=0. (4) 

These relations correspond to those in Ising model: 

[ot,<Tj<Tj + i] = 2afajaj +1 (5i !j +5 l j+i), 
[of,oj] = 0, 

[« +1 ,a^| +1 ]=0. (5) 



Then we obtain the mapping 

Fi^af, t? afaf +1 . (6) 

Accordingly, the Hamiltonian (2) can be mapped onto 
the ID quantum Ising model like following: 

H 'w -> - zZ XX^z + J<i<?ai+i), (7) 

ol i 

where the subscript a implies there is one or more Ising 
chains. The number of Ising chains is determined by the 
size of the square lattice for the Wen-plaquette model 
(see Ref. 11 ). Since the Ising chains decouple from each 
other, we can consider only one Ising chain without loss 
of generality, and reduce the Hamiltonian (7) to 

H w = -JY&wot + 9w = j. (8) 

i 

Then we may explore the quantum properties of the 
original Wen-plaquette model by studying the corre- 
sponding ID Ising model. 



III. STRING ORDER PARAMETERS IN 
WEN-PLAQUETTE MODEL 

For the ID transverse Ising model (3), there are two 
phases: in the limit of gi 3> 1, the ground state is a para- 
magnet with all spins polarized along x-axis, (erf) — > 1; 
in the limit of gi <C 1, there are two degenerate ferro- 
magnetic ground states with all spins along positive or 
negative z-axis and (erf) — > 0. Consequently, in Hamilto- 
nian (8), there is a quantum critical point at 21 

gw = j = l (9) 

that divides the two phases. Accordingly, the original 
transverse Wen-plaquette model also has two phases sep- 
arated by this quantum critical point — in the region of 
gw > 1; the system is a topologically ordered state; in 
the region of gw < 1, it's a spin-polarized state. 

Noting that the local order parameters cannot be used 
to learn the nature of TQPT any more, we introduce 
two non-local order parameters "01 and 1^2 as string or- 
der parameters (SOP's) in the transverse Wen-plaquette 
model. They are defined by the expectations of string op- 
erators ]J i Fi and ]J i rf with i as the site index along a 
string in the diagonal direction 11 , i.e. ip\ = (Hi ^i) an d 
i'2 = (Hi-rf), respectively. 

We then calculate these two SOP's by using the map- 
ping in Eq.(6). For -02, one has 

^ = (ll T ^ = (°i<>Z<4<4 ■ ■■<-!<) = («) , 

(10) 
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which becomes the correlation of two spins in the Ising 
chain of length n. By employing the Jordan- Wigncr 
transformation 

< = 1 - 2 4c«, < = -(c + (11) 



where 



!/ = n (i-24cm) = n ( c 



where .4 



4 



B — rt 



?71 <-m 5 ^ m 



(12) 

, cj„ and c rn are 



the creation and annihilation operators for fermions, we 
obtain 

V> 2 = K<7*) = <J(ci + cl)i/(c + )) 

= (S X i4 2 S a • ■ ■ B n _ x A n ) . (13) 

Following the Wick's theorem, we can transform Eq.(13) 
into a Toeplitz determinant as 22 



G\2 G13 • • • G\n 

G22 G23 ■ ■ ■ Gin 

G n -12 G ti -13 • • • G n -\ n 



(14) 



where 



° 11 "27 



" dh g w -cosk-isink plk{t _ j} 
* J (gw - cos k) 2 + sin 2 k 



(15) 

In the thermodynamic limit N — > oo, we have (see Rcf. 23 
or Ref. 24 ) 



Ij'2 



(! - 9w) 1/4 when 9w < 1 
when c/vk ^ 1 



(16) 



The exponent 1/4 here agrees with the critical exponent 
2f3/v = 1/4 (see Ref. 21 ). 

By the same method we can obtain the result of if)\, 



IV. QUENCH DYNAMICS OF TQPT 

In this section we study the dynamics of TQPT in 
the transverse Wen-plaquettc model. The Kibblc-Zurck 
mechanism (KZM) 25,26 is a general theory to explore 
the dynamics of second order phase transitions including 
the quantum case. According to KZM, during a quench- 
induced phase transition, the system undergoes three 
stages of evolution: adiabatic-pulse-adiabatic. It pre- 
dicts that the density of topological defects, which are 
generated by the pulse evolution, is a function of quench 
time. 

We can solve the dynamic problems in the Wen- 
plaquette model by taking a sequence of transformations 
as in Ref. 27 . First, through the Jordan- Wigner trans- 
formation (11), the spin operators a*' z are represented 
by fermionic operators c„. Second, the operators c„ are 
Fourier transformed into the momentum space: 



k 



c k e 



(20) 



(21) 



Third, after Bogoliubov transformation, one have 
Cfe = u k r] k + v*_ k rji k , 

where rjk and ry k are fermionic operators. For dynamic 
problems, the expression (21) becomes 



Cfc(i) = u k (t)rj k + vt k (t)r]l k . 



(22) 



Now, a dynamic solution can be written in the form of 
Bogoliubov mode (u k (t), v k (t)) after this whole transfor- 
mation procedure. 

By this method, the dynamics of the quantum Ising 
model can be expressed as the time evolution of the Bo- 
goliubov mode via so-called Bogoliubov-de Gennes dy- 
namic equations (see Ref. 27 ): 



(17) 



ihdu k /dt = +2(gi(t) — cos k)u k + 2 sin kv k 
ihdv k /dt = —2(gi(t) — cosk)v k + 2smku k 



(23) 



In addition, we can take advantage of the duality of 
ID Ising model (see Ref. 22 ) : 



<of +1 and sf = Y[ <4 



(18) 



k<i 



Thus the result in Eq.(17) is turned into tpi = to, where 
m = (erf) is the spontaneous magnetization. We can find 
that when N — > oo. the spin correlation (10) is the square 
of m. Consequently, from Eq.(16), we have 



0! 



when gw < 1 

(l-.9w) 1/8 when g w > 1 



(19) 



From the calculations above, one may notice that the 
non-local SOPs in a 2D transverse Wen-plaquette model 
are transformed to the local order parameters in the dual 
ID Ising model. 



Consider a linear quench, i.e. 



0) = , 



(24) 



where t varies from — oo to and the quench time tq 
characterizes the quench rate (defined by 1/tq). Equa- 
tions (23) can be transformed into the form of Landau- 
Zener (LZ) model 28 (the connection between the KZM 
and the LZ model can be found in Ref. 29 ' 30 ) : 



ihdu k /dr = 
ihdv k /dr = 



-i(rA fc )ufc + i«fe 
-^(rA k )v k + jiifc 



(25) 



where 



r = 4tq sin k( h cos k), 
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FIG. 2: (Color online) The expectation value of Fi at 
t = varying with the quench time tq, where /i = 1 — 
IJ_* i( lfc(i±^-c- 2 "« 

h = —m. 



e 27TT Q s,n fc ) and the approximative result 



and 



A fc 1 = 4tq sin 2 k. 



From equations (25), we derive a second order differ- 
ential equation for v k : 



d 2 v k 



dr 2 

After the substitutions 
1 



iAi- 1 . 
+ (-T 2 A 2 k + ^ + -H = 0. 



4iA* 



and z 



t/4 



we have 



1 



z 2 )v k (z) = 0. 



(26) 



(27) 



(28) 



This kind of differential equation has a general solution 

Vk(r) = -[aD- s - 1 (-iz) + 6D_ a _i(iz)], 



Ufc(r) = i-A k T + 2i£)v k (T), 



(29) 



where D m (x) is the so-called parabolic cylinder function 
(PCF) or Weber- Hermite function 31 . 

According to the boundary conditions and the char- 
acters of the PCF, one may derive the approximative 
solutions to equations (28) at the end of linear quench 
fori^O (SeeRef. 27 ): 



1— cos k 



_|_ g Itttq sin 2 k 



|Ufc| 

\v k \ 2 = 1 - \u k \ 2 , 
u k v k = I sm ^ + sgn(k)e 



tttq sin 



g— 7TTQ sin 2 k^ifk 

(30) 

with the condition tq > 1. 

Now we calculate the dynamic solutions in the Wen- 
plaquette model via the same transformation procedure 
plus the mapping (6). We shall also consider 

g w (t<0) = - — . (31) 



The quenching process can be set as tuning the strength 
of the transverse field from zero to very large compared 
with the coupling constant g. The phase transition is thus 
from a topologically ordered state to a spin-polarized 
state during the time-evolution from t — > — oo to t = 0. 

First, we calculate the expectation value of Fi after the 
quenching process. From the Jordan- Wigner transforma- 
tion, er* = 1 — 2c\ l c n , we have 



(Fi. 



<(l-2ctc„)) 



Through the Fourier transformation, we express & m c n in 
momentum space 



„t 



1 

27 



dk 



dk'clc k ,e i{k ' n ~ km) . 



(32) 



After the Bogoliubov transformation (21), we obtain 



:!2 



— / dfc 

2tt 



1 

27 



T 

dk\v k \ 2 e lk{n - m) . (33) 



From the solutions (30), we derive 



(Fi, 



1 - 2 

— Itzt q sin 2 k\ 



l 

27 

TQ»1 



dfc( 



1 + cos k 



(34) 



We can see the original state in the limit of g ^> 1, 
which is a topological ground state with (Fi) — > (a*) 
1, will finally evolve into the trivial spin-polarized state 
with (Fj) — > (a*) — > 0. Figure 2 shows how the expecta- 
tion value of Fi at the end of quenching process (t = 0) 
varies with the quench rate. 

We can obtain the the numbers of Zi charges and Z2 
vortices by defining 



and 



m d = \ (1-^)= E c l< 



(35) 



zGodd 



i£odd 



respectively. Before the quench, the system is in a topo- 
logical ground state. There is no quasiparticlcs (N c = 
(jV c ) = 0, N v = (Af v ) — 0), so wc have Fj = 1 for all pla- 
quettes. At the end of the quench, we can also calculate 
the density of plaquettes Fi — — l 33 by 



dk( 



1 + cos k 



27TTQ sin 2 k\ TQ ^ 1 I L 



2n,/2T < 



Q 



(36) 



FIG. 3: (Color online) The colored curves represent |«fe(r)| 2 in 
the original form of the PCF varying with momenta k from — tt 
to 7r, where the red one is at t = — 0.5tq, the blue at t = —tq, 
the green at t = —2tq with tq = 50 in all the three cases. The 
dashed, dotted and dot-dashed curves, respectively, represent 
the approximate function for |vfc(r)| with corresponding pa- 
rameters, with the negative parts under the x-axis being cut 
off. They match with the colored curves very well 



It is obvious that when tq — > oo, half of the plaquettes 
will be turned into Fi = — 1, which implies (Fi) — > 0. 

However, this result cannot give us any information 
about the quenching process or the critical behaviors. 
In order to obtain such information, we find that, for 
tq 2> 1, the expression 



M*)i 2 = iki+ 



cos k + t/rq 



^/l + 2t/TCOSk+{t/T Q ) 



f)- 



-2-tttq sin k 



(37) 

is a time-dependent approximate function for the general 
solution |tifc(r)| 2 in Eq.(29), if we cut off the negative part 
of the curve (see Fig. 3). 

By using this approximate function, we calculate the 
value (Fi) during the quenching process as a function 
of the time t near the quantum critical point. The re- 
sult is shown in Fig. 4. Further calculations show that 
the derivative of the expectation- value (Fi) diverges at 



point t = t c 



-TQ 



dt 



-co, which characterizes 



the critical point. In particular, we obtain a logarithmic 
scaling law of quench dynamics near the quantum phase 
transition as 



d(F t ] 
dt 



]n\t-t c \. 



(38) 



Such a dynamics is analogous to the static scaling behav- 
ior of the specific heat near the critical point for the ID 
quantum Ising model. 

The physical picture can be described as follows. At 
beginning, the external field is weak and can be treated 
as a perturbation, which creates vortices or charges and 
drives them to move around and annihilate each other. 




FIG. 4: (Color online) The expectation value of Fi during the 
quenching process varying from t = —2tq to t = with tq = 
50 and the absolute value of the slope (| d |) corresponding 

to the curve above. The small graph shows | d ffi | ~ In |i/ T Q + 
1| from both sides approaching the critical point t = —tq. 



As a result, the density of quasiparticles remains small. 
As the strength of the field continuous to grow, the den- 
sity of plaquettes Fi = — 1 increases rapidly near the 
point gw = 1 (that is t = — tq). The rate of density- 
changing diverges with a logarithmic scaling law. Finally, 
at the end of the quench, the field becomes so strong that 
the vortices and charges are all confined and cannot be 
treated as quasiparticles. If the quenching process is in- 
finitely slow, half of the plaquttes overturns. 

In addition we shall notice that only the linear quench 
is considered here. However, for other cases, our method 
is not reliable. Inspired by other papers on quantum 
quench in the toric-code model (e.g. Rahmani et a/. 34 , 
where a sudden quench is studied), we can study the 
quench problems more generally through the time evolu- 
tion of the entanglement entropy in the Wcn-plaqucttc 
model. Nevertheless, we will not consider this in this 
work. 
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V. CONCLUSION 

In summary, we study the dynamics of TQPT caused 
by a linear quench in the transverse Wcn-plaqucttc 
model. We first show how to derive the mapping from the 
2D Wen-plaquette onto the ID Ising model by compar- 
ing their commutation relations. Based on this mapping, 
we point out the quantum critical point in the transverse 
Wen-plaquette model and calculate its non-local order 
parameters. We then calculate the expectation value of 
Fi at the end of the quench, and further show how this 
value varies during the whole process. In particular, we 
find a logarithmic scaling law of quenching process near 
the TQPT. 

Finally we address the realization of the Wen-plaquette 
model in an optical lattice of cold atoms. Because the 
Wcn-plaqucttc model can be regarded as an effective 
model of the Kitaev model on a two dimensional hexag- 
onal lattice, one may first realize the Kitaev model. The 
Hamiltonian of the Kitaev model is 18 



The key point is to calculate [Fi,Fj], 

= W X T V T X - ~T V T X T V T X - ~T V 1 
L'l i+x i+y+x i+y ' 'j j+x j+y+x' j+yl 

= W X T X T V T X - ~T V ]t V T X - ~T V +T X 

L 1 ' J j+x J+y+x j+yl i+x ' i+y+x i+y 1 '? 

\ T v T x y x _ y i x _ y , x y 
l i+x' j j+x J+y+x j+yl i+y+x' i+y 1 'i 'i+x 

W X - - T X T V T X - ~T y \T V +T X T V T X - - 

I i+9+i' J j+x 1 j+y+x 1 j+yl' i+y ' i i+x ' i+y+x 

\ T V T x y x y I ( 41 n 

L'i+i' 'j ' J+x 1 J+y+x 'j+yl- \^ L ) 

According to the commutation relations of Pauli opera- 
tors, 



W X T X T V - 
l'l ' 'j 1 J+x 



j j+X j+y+X j+yl 

C T X \t V T X - ~T 
' J 1 j+x J+y+x j+y 

t-x y r x x i y , 
T j T j+xVi ' T j+y+i\ T j+y "t 

■ ZIO^j+xTj T t T J+ y +s .T j+ . 



T X W X T y \T X - ~T V 

J L't ' ' j+il' J+y+x' j+y 

, ~W X T V 1 
-y+xl't ' 'j+yl 

2i() -i ~T X T y T X - ~T Z - 
■""i.J+y'j 'j+i 'j+y+x 'i ' 

(42) 



T X T V -T X 
'j 'j+x'j- 



H = ^2 i^ a j,l a j+i,i 

j-\-l=even 



(39) 

where j and I denote the column and row indices of the 
lattice. In the limit of 3 Z ~ 3 y in this model, the 
effective Hamiltonian of Kitaev model is simplified into 
that of the Wen-plaquette model as 



T V X V x y 1 

'i+x' '3 j+i j+y+x j+yl 

= \t V T X W V T X - -T V 
I' i+x' J 1 j+x ' j+y+x j+y 



T X T V -\t V ~ T 
' 1 ' i+9. Y'i.+t.i ' ■ 



X 



'j ' j+x L ' i+x ' ' j+y+il T j+y 



u 



-9? A- , - T Z T V T X - -T° 



xr y y i x y 
'j I' i+x' 'j+il' j+y+x' j+y 

T X T V T X - ~W V T V 1 

J j+x ' j+y+xl' i+x' ' j+yl 

UK- ■ , ~T X T V T Z . - 

^ lu i,3+y'j ' j+x' j+y+x 



j+y' 



(43) 



H ° ~ ~ 16|J I 3 ^ a i clt ( l ) a * i shm (J upW (J l 



own(i) ' 



(40) 



Then one can use the Kitaev model in the limit Jx^> Jz ~ 
J y on a torus to do the TQC. The realization of the Ki- 
taev model on the 2D hexagonal lattice has been pro- 
posed in Ref. 35,36 . The essential idea realizing the Kitaev 
model is to induce and control virtual spin-dependent 
tunneling between neighboring atoms in the lattice that 
results in a controllable Heisenbcrg exchange interaction. 
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Appendix 

In this part, we give detailed calculations about com- 
mutation relations (4), which is related to the consistency 
between Wen-plaquette model and Ising model. 



r x T v t x - ~T y 1 

3 J+x' J+y + X 'j+yl 

- W X - - T X W V T X - -T V +T X \T X . . T V W X - - 
L't+x+x) 'j 1' j+x' J+y+x ' j+y ' 'j l'i+x+x' ' j+il' j+y+x 

T V +T X T V W X - - T X - .W V +T X T y T X - - 
J+y T 3 j+xV l+x+x' 'j+y+xi'j+y ^ ' J > j+x' J+y+x 



W X - - T V 1 

L' i+x+x' j+yl 



2iSi 



+V,j T j T j+x T j+y+x T j+y 1" ^ l °i+i,j T j ^j+x T j+y+x 



'j+y 



(44) 



T V X V x y 
i+y' J j+x' j+y+x' j+y 

-V T x^V ^x ^y 
i+y' 3 

rV n-X 1J 



J i f,' T i \ T j+i T 3+y+i T j+y + T j Vi+y' T j+il T J+y+x T 



T j T j+x l T i+y' T j+y+i\ T j+y 



x\ T V T V } T x . _» 
3 t'l+y' j+il' J+y+x' j+y 

y n-y i 



Hi) , - T Z T V T X . -T V 
^ L °i+y,J T j T 3 +x T J+y+x T j+y 



T J T 3+x T 3+y+xi T i+y' T j+y 

v X ~-V 



^ iu l,3+xlj 'j+x' j+y+x 
Tf-H' ( 45 ) 

we can calculate (41) in following four cases: 
1. when i = j + x, 

[Fi,Fj] 



- 7iT X T Z -T X - ~T V T V T X ~ ~T V 

- &t,l 3 I J+x' J+y+X 1 j+y 'i+X 1 t+V+x' i+y 

T X T V T Z ~ ~T V 
J j+i J+y+x j+y 

= -1T X T Z -T Z - -T V T V T X - 

3 J+x' j+y+x' j+y' i+x %+y+x 
z y 
T j+y+x 7 j+y 



oj T x y x 

'i+i'i+y+x 



' i+i' i+y+x j 'j+x 



0; 



(46) 
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It is clear that we also have [Fj, Fj] — in other cases. 

2. when i = j + y, 

3 T j+x T j+y+x T j+y T i+x T i+y+x T i+y ~ T j T j+x Then, we only need to verify [i^,rj] = 2FiT x (Sij- X + 



~ LU ] 'j+x 'j+y+x 'j+y 'i+x ' i+y+x 'i+y zt 'i 'j '> 

t z , - , -t v ~t x , - , -r y ~ 5i i-4j), while others can be easily obtained from the com 

' 3+y+x' j+y' i+y+x' i+y l <-> yn J 



y 

T i+y+x T i+y 



T X T V T Z ~ ~T Z ~T X ~ ~T V + 2"T? ~T X T V T Z - - 

'j j+x j+y+x 3 +y %+y+x i+y ' ^'j+y'j 'j+x'j+y+x shown in the following 



mutation relations of Pauli operators. The verification is 



(47) 



3. when i + x = j 



-7iT X T Z T V T X - ~T V T X - ~T V +2iT X T V T X T V 
i 3 j+x 'j+y+x 'j+y 'i+y+x 'i+y ~ lLl '% 'i+x' 3 'j+x 



[Fi , Tj } 

3+y+x i+y i+y 

- o T x z y x z y ,r,xz y x = ~ ^10 j^i+x^i ^i+ x T i+y+x T i+ v ~ ll6 j^+y' i ' i+ x ' i+y+x' i+y 



T Z ~T V 

' J+y 'i+y 

= 0; (48) 
4. when i + y = j 

— 2-!T X T V T X T Z -T X - ~T V T V -2lT X T V T X - - T Z 
i4/ J 'i+x'j ']+x']+y+x'j+y'i+y i+x i+y+x j 

T V T X - -T V 

j+x i+y+x 'i+y 

= -2T X Tl A T Z T Z T X ^, + 2T X Tl A T Z 



;T; T- 



i 'i+x' j j+x j+y+x j+y 1 *"* ' i+x' i+y+x ' j ' j+y+x 

V 



T i+y 



(49) 



^3,1+xTi T i+s .T i+ -T l+y+i .T i+ - Z,0 J}l +yT l T i+x T l+ y +x 



T X ~T V 

't+y 'i+y 



OA t x P 9A t x P 

^ u j,i+x'j r i j^+y j i 

2F i if(6 jti + S! + 6j !i +$). (50) 
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